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Abstract 

The non-equilibrium self-consistent generalized Langevin equation theory of colloid dynamics 
is used to describe the non-stationary aging processes occurring in a suddenly quenched model 
colloidal liquid with hard-sphere plus short-ranged attractive interactions, whose static structure 
factor and van Hove function evolve irreversibly from the initial conditions before the quench to a 
final, dynamically arrested state. The comparison of our numerical results with available simulation 
data are highly encouraging. 
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I. INTRODUCTION. 



The non-stationary, slowly-evolving dynamics of deeply quenched fluids, referred to as 
aging, has been the subject of considerable attention over the last decade [l[ Concen- 
trated emulsions {3], colloidal gels {4], and aqueous clay suspensions J5] are some examples of 
aging systems. In spite of the apparent diversity of these structurally-disordered and out-of- 
equilibrium materials, the appearance of certain universal features in their non-equilibrium 
evolution suggests the existence of an underlying common source of the observed dynamic 
properties. Although this non-stationary behavior is associated with the formation of dis- 
ordered solids, including hard materials such as polymer glasses p, the main features are 
best exhibited by soft materials such as those above. In particular, the study of the dynamic 
properties of aging colloidal glasses and gels is specially interesting, since the observations 
provided, for example, by experimental methods such as dynamic light scattering [t lo| can 
sometimes be complemented with direct visualizations at the level of individual particles 
by means of digital video imaging techniques ll|, [l2| . Computer simulation experiments in 



well-defined model systems have also contributed with important complementary informa- 
tion about the general properties of aging [13l4l5]. 



From the theoretical side the study of aging has been addressed in the field of spin 



glasses, where a mean-field theory has been developed within the last two decades [16]. 
The models involved, however, lack a geometric structure and hence cannot describe the 
spatial evolution of real colloidal glass formers. About a decade ago Latz Jl7] attempted 
to extend the mode coupling theory (MCT) of the ideal glass transition 18|, |l9| to describe 
the irreversible relaxation, including aging processes, of a suddenly quenched glass forming 
system. Similarly, De Gregorio et al. [20| discussed time-translational invariance and the 
fluctuation-dissipation theorem in the context of the description of slow dynamics in system 
out of equilibrium but close to dynamical arrest. Unfortunately, in neither of these two 
theoretical efforts quantitative predictions were presented that could be contrasted with 
experimental or simulated results in specific model systems of structural glass-formers. 

For concreteness, let us focus our discussion on the conceptually simplest glass-forming 
system, namely, a mono-component fluid made of N identical spherical particles in a vol- 
ume V which interact through the pair potential u(r) (although in the experimental realiza- 
tion of this idealized model we probably have to consider a small amount of polydispersity 
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to suppress the kinetic pathway to ordered phases). Assume that in the absence of ex- 
ternal fields this system is initially prepared in an equilibrium state corresponding to a 
mean density = N/V and a temperature T^°\ in which the static structure factor is 
S^°\k) = S eq (k;n,T^). In the simplest idealized quench experiment, at the time t = 
the temperature of the system is instantaneously and discontinuously changed to a value 
Ttf>. Let us assume that along the process that follows the quench, the density and the 
temperature are constrained to remain uniform and constant, i.e., that n(r, t) = and 
T(r, t) = T^' at any position r in the volume V and any time t > 0. The relevant question 
then refers to the value of the static structure factor S(k; t) for t > 0, and to the evolution 
of the dynamic properties of the system along this process. 

The referred dynamic properties can be described in terms of the relaxation of the fluctu- 
ations Sn(r, t) of the local concentration n(r, t) of colloidal particles around its bulk equilib- 
rium value n. The average decay of 5n(r, t) is described by the two-time correlation function 



F(k,r;t) = iV 1 5n(k, t + r)5n{— k, t) of the Fourier transform 5n(k, t) of the fluctuations 



5n(r, t), whose equal-time limit is S(k;t) = F(k,r = 0;t) = N~ 1 5n(k, t)5n(— k, t). We refer 
to the time r as the correlation time, and the overline refers to the average over the prob- 
ability distribution of the non- equilibrium ensemble that governs the statistical properties 
of d~n(r, t) at the evolution time t. This ensemble will surely coincide with an equilibrium 
ensemble only in the limit t — > oo, provided that no dynamic arrest condition appears along 
the process. 

After the sudden temperature change at t = has occurred the system evolves sponta- 
neously, searching for its new thermodynamic equilibrium state, at which the static structure 
factor should be S eg (k;n^°\T^). If the end state, however, is a dynamically-arrested state 
(a glass or a gel), the system may never be able to reach this equilibrium state within experi- 
mental times; one then refers to the evolution time t as the waiting or aging time [l, 8, 10-12]. 



The dependence of S(k; t) and F(k, r; t) on t characterizes the non-equilibrium evolution of 
the system, whose quantitative theoretical first-principles description, to our knowledge, has 
not been available until now, in spite of important theoretical efforts like those referred to 
above. 

In recent related work [2l|, however, an extension was proposed of the self-consistent 
generalized Langevin equation (SCGLE) theory of colloid dynamics 22M26I] and dynamic 



arrest 



27H34]. aimed precisely at describing this non-equilibrium evolution of S(k;t) and 
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F(k,r;t). This extension was based on Onsager's theory of thermal fluctuations 
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adequately extended 



40 



41] to allow for the description of memory effects. The purpose of 



the present paper is to provide the first practical and concrete application of such general 
non-equilibrium theory of colloid dynamics by means of its use in the quantitative description 
of the aging process of a model monocomponent glass-forming liquid. 

In this particular context, such non- equilibrium self-consistent generalized Langevin equa- 
tion (NE-SCGLE) theory consists of a closed self-consistent system of equations for S(k; t) 
and F(k, r; t), which we numerically solve here for a model mono-component fluid of parti- 
cles interacting through the hard-sphere plus short-ranged attractive Yukawa potential. This 
model system exhibits the glass-fluid-glass reentrance predicted by the equilibrium SCGLE 
theory [30] (and originally discovered by MOT 42|). Here we discuss the isochoric quench 
of this fluid from an initial equilibrium state (n^°\T^) in the reentrant fluid pocket of the 
(n, T) state space, to a final temperature T^> (< T^) in the vicinity and below the attrac- 
tive glass transition temperature r"(n' ') corresponding to the density n^°\ This process 
mimics the computer simulation aging experiment reported by Foffi et al. in a similar 
model system (hard-sphere plus short-ranged square well). Here we discuss our theoretical 
predictions in reference to the observed behavior in this simulated quench experiment. 

We start this discussion by summarizing in the following section the full non-equilibrium 
self-consistent generalized Langevin equation (NE-SCGLE) theory, which does not involve 
the restrictive assumption of spatial homogeneity. In the same section this theory is simpli- 
fied according to the assumption that the system is constrained to remain spatially homo- 
geneous and isotropic. The actual solution of the resulting equations are reported in Sec. 
IIHI The last section contains a summary of the results. 



II. NON-EQUILIBRIUM SELF-CONSISTENT GENERALIZED LANGEVIN 
EQUATION THEORY 

The previous discussion implicitly assumes that S(k; t) and F(k, r; t) adequately represent 
the structural and dynamic properties of the quenched system along the irreversible equili- 
bration process. This assumption, which is thought to be accurate in the absence of external 
fields, or when the effects of these external fields are very small, is in reality a strong simpli- 
fying assumption when the local concentration fluctuations do not relax within experimental 
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times, as it occurs at and near dynamically arrested states. The general non-equilibrium 



SCGLE theory proposed in ref. 



21], however, does not incorporate this simplifying assump- 



tion at the outset. Instead, it describes the statistical properties of the instantaneous local 
concentration profile n(r,t) of the colloidal liquid in terms of the coupled time evolution 



equations for its mean value n(r, t) and for the covariance cr(r, r'; t) = 5n(r, t)5n(r', t) of the 
fluctuations 5n(r,t) = n(r,t) —n(r,t). In this section we briefly review the general NE- 
SCGLE theory and then particularize it to instantaneous homogeneous quench processes. 



A. General NE-SCGLE theory. 



The referred equations for n(r, t) and a(r, r'; t) read 21] 

dn{r, t) 



()f - D°V • 6(r, t)n(r, t)V^[r; n(t)} (2.1) 
and 

daij^yt) =jD o v . - (r; t) fo(r ^ )y I driS[rjri;n (t)} a ( ri ,T';t) 

+ D°V • n(r\ t) b{r', t)V J dr 1 S[r', r i; n(t)]a(r 1 , r; t) ( 2 - 2 ) 

- 2D°V ■ n(r, t) b(r, t)V5(r - r'), 

in which Dq is the self-diffusion coefficient of the colloidal particles in the absence of direct 
interactions, /x[r; n] is the electrochemical potential at position r (which is a functional of 
the local concentration profile n(r)), and £ [r, r'; n(t)] is the functional derivative £[r, r'; n] = 
[5(3fi[r;n]/5n(r')] evaluated at the concentration profile n(r) = n(r, t). Thus, for given Dq 
and fi[r; n], Eqs. (12. lft and (12. 2 p would constitute a closed system of equations for n(r, t) and 
<r(r, r'; t) if it were not for the presence of the dimensionless local mobility function b(r, t). 

This mobility function b(r, t) describes the local frictional effects of the direct (i.e., conser- 
vative) interactions between the colloidal particles, as deviations from the value b(r,t) = 1, 
and can be expressed in terms of the memory function of the two-time correlation function 



C (2) (r,r';t,t') = 5n(r, t)8n(r', £'), which we write as C (2) (r, r + x; t, t + r) = C(x,r;r,t), 
and which is the value of the van Hove function at a spatial location r in the system and at 
an evolution time t. Since the covariance is <r(r, r'; t) = C^(r, r'; t, t) = C(x, r = 0; r, t), it 
can also be written as cr(x; r, t). Although in the development of the non-equilibrium SC- 
GLE theory the assumption of absolute spatial homogeneity and isotropy is avoided, these 



spatially-varying van Hove function and covariance do depend on the location r in space 
but are assumed to be approximately isotropic within a small volume around r, so that they 
only depend on the magnitude |x| of the correlation vector x. Under these conditions, the 
local covariance c(|x|; r, t) can be written in terms of its Fourier transform a(k; r, t) as 

a{\ x |;r,t) = J d 3 k e -^a(k;r,t), (2.3) 

so that Eq. (12. 2p may be re-written as 

9a( ^ r,t) = - 2k 2 D°n(r, t)b(r, t)£(k; n(r, t))a(k; r, t) 

(2.4) 

+ 2k 2 D°n(r,t) b(r,t), 

where 8(k; n(r, t)) = (27r)~ 3 J <i 3 /ce _ik ' x £[r, r + x; n(r, £)]. Similarly, the local van Hove 
function C(|x|, r; r, t) can also be expressed in terms of its spatial Fourier transform as 

C(|x|,r;r,t) = -i- J d 3 ke-^C(k, r; r, t). (2.5) 

Let us notice that we can also introduce the notation C(k, r; r, t) = rT(r, t)F(k, r; r, t), with 
F(k,r;r,t) being the non-equilibrium intermediate scattering function, whose initial value 
F(k, t = 0; r, t) = S(k; r, t) defines the time-evolving spatially-varying static structure factor 
S(k; r, t); this more familiar notation will be employed later on. 



According to Ref. (21] . the actual calculation of the local mobility function 6(r, t) requires 
the solution, at each position r and each evolution time t, of a system of equations involving 
the Laplace transform (LT) of C(k, r; r, t) (denoted by C(k, z\ r, t) = J °° dre~ ZT C(k, r; r, t)), 
as well as the LT of its self component Cs(k, r; r, t), and of the r-dependent friction function 
r,t), namely, 

C{k, Z] r, t) = k *D°n{r,t)<!-\k;r,t) i ( 2 ' 6 ) 

Z + l+\{k) A(*(z;r,t) 



C s {k, z; r, t) = ^5 , (2.7) 

z H — 

1+A(fe) A(*(z;r,t) 

and 



AC*(r;r,t) = D ° I dk k 



3(2 



7T 



a(k; r, t)/n(r, t) — 1 
a(k;r,t) 



2 



C(k,T-r,t)C s (k,T-r,t). (2.8) 
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with X(k) being a phenomenological "interpolating function" given by 

X(k) = — 1 



21 



28 



22| 



1 



2 ' 



(2.9) 



where k c > 2n/d, with d being some form of distance of closest approach. A simple empirical 
prescription is to choose k c as k c = k min , the position of the first minimum (beyond the main 
peak) of the static structure factor S(k;r, t) = a(k;r,t)/n(r,t). The local mobility b(r, t) 
finally follows from the solution of these equations by means of its relation with A(*(z; r, t), 
namely, 



b{r,t) = l + AC*(z = 0;r,t) 



n -1 



(2.10) 



B. Instantaneous homogeneous quench. 

Let us now discuss the application of this general theory to the particular conditions 
referring to the irreversible evolution of the structure and dynamics of a system constrained 
to suffer a programmed process of homogeneous compression or expansion (and/or of cooling 
or heating). Under these conditions, rather than solving Eq. fl2.ll) for n(r;t), we assume 
that the system is constrained to remain spatially uniform, n(r;t) = n(t), according to a 
prescribed time-dependence n(t) of the uniform bulk concentration (and/or to a prescribed 
uniform time-dependent temperature T(t)). In consistency with this assumed constraint we 
have that the dependence on the position r disappears from the previous equations so that, 
for example, Eq. (12 .4p may be re-written as 

d(T ^ Q = -2k 2 D°n(t)b(t)S{k; t) [a{k; t) - £-\k- 1)] , (2.11) 

with S(k; t) = S(k; n(t)) and with 



bit) 



drAC(r;t) 



o 



(2.12) 



where A£*(r;t) is provided by the solution of the self-consistent system in Eqs. f l2.6l) - fl2.8p 
for the uniform bulk concentration n(r; t) = n{t). 

Among the many possible programmed protocols (n(t), T(t)) that one could devise to 
drive or to prepare the system, the simplest corresponds to the idealized quasi-static process, 
in which the relaxation rate da(k; t) / dt is virtually negligible due to a virtually instantaneous 
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"thermalization" of a(k;t) to its local equilibrium value a Le '(k;t) = l/£(k;t) [2l|, |39|. A 
quasi-static process, however, is a rather unrealistic concept, at least in the limit of small 
wave- vectors, in which the relaxation times diverge as k~ 2 , as seen in the example below. In 
contrast, a far more interesting and fundamental protocol corresponds to the opposite limit, 
in which the system, initially at an equilibrium state determined by initial values of the 
control parameters, (n^°\ T^), must adjust itself in response to a sudden and instantaneous 
change of these control parameters to new values (n^,T^), according to the "program" 
n(t) = n^9(-t) + n^9(t) and T(t) = T^9(-t) + T^6(t), with 6(t) being Heavyside's step 
function. 

Under these conditions the formal solution of Eq. (12. lip can be written, for t > 0, as 



a 



(k;t) = a°(k)e- a{k)u{t) + [£ ^{k)]" 1 (l - e - a{k)u{t) ) 



(2.13) 



where £^'(k) = £(k;n^' ,T^') is the Fourier transform of £[\ r — r' \;n^' ,T^'] 
[5Pfi[r;n}/5n(r')} n= - (f)tT=T(f)j 



u(t) = / b(t')dt', 



(2.14) 



and 



a(k) = 2k 2 D°n {f) £ if) {k). (2.15) 

Clearly, the presence of the time- dependent mobility b(t) couples this formal solution with 
the self-consistent system in Eqs. (I2.6l) - (l2.8p . For the present conditions, and in terms of 
the non-stationary static structure factor S(k;t) = a(k;t)/n^ and intermediate scattering 
function F(k, r; t) = C(k, r; t)/n^\ we may rewrite such self-consistent system of equations 

as 



F(k,z;t) 



S(k;t) 



k 2 D°S- 1 (k;t) ' 



Z + 



(2.16) 



1+A(fc) A(*(z;t) 



F s {k,z;t) 



z + 



k 2 D ) 



(2.17) 



1+A(fc) A(*(z;t) 



and 



AC(r;t) 



Dn 



dk k 



S{k;t) - 1 
S{k;t) 



F(k,T-t)F s (k,r;t). 



(2.18) 
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Eqs. fl2.12p - ()2.18p constitute our general self-consistent description of the spontaneous evo- 
lution of the structure and dynamics of an instantaneously and homogeneously quenched 
liquid. 

Of course, one important aspect of this analysis refers to the possibility that the end state 
of the quench process happens to be in the region of dynamically arrested states. For the 
discussion of this important aspect it is useful to consider the long-r (or small z) asymptotic 
stationary solutions of Eqs. fl2.16p - fl2.18l) above. Just like in the equilibrium SCGLE theory 



281 ]. these may be analyzed in terms of the asymptotic values of these dynamic properties 



(the so-called non-ergodicity parameters), given by [21] 



fit t) = lim F ^ T]t) = (2 19) 

- S(k) X(k;t)S(k;t) + k^(t) 1 ^ 



T— >OQ 



and 



MM S & F «W = W+%«ty (2 - 20) 



where the squared localization length 7(2) is the solution of 



1 r dkk* J^l:^3 M ! „,„. (2.21) 



7(t) 6tt 2 ^ J [\(k; t)S(k; t) + k 2 j{t)} [A(fc; t) + k 2 -y(t)] 

These equations are the non-equilibrium extension of the corresponding results of the equilib- 
rium SCGLE theory, and their derivation from Eqs. f l2.16p - fl2.18p follows the same arguments 



as in the equilibrium case [24 1. 



In the following section we numerically solve Eqs. fl2.12p - fl2.18p for still more specific 
conditions, namely, for an isochoric quench of a model colloidal system, in which n^' = 
and Eq. (I2.13P can be written in terms of the time-evolving static structure factor S(k; t) = 
a(k;t)/n( f) as 

S(k; t) = S\k)e- a{k)u{t) + Sf{k) (1 - e - a ( fc M')) , (2.22) 

with Sj 9 (k) = [n^'£«)(fc)] -1 . Let us notice that in the limit in which the friction function 
A(*(r;t) vanishes, b(t) = 1 and hence u(t) = t, so that Eq. f)2.22p reads 

S*(k; t) = S°(k)e- a{k)t + S e /(k) (1 - e~ a{k)t ) . (2.23) 

This limiting expression describes an exponential interpolation of S(k; t) between its initial 
value S°(k) and its final equilibrium value Sj q (k) = [n^5^(fe)] _1 . It is then important 
to notice that the general solution S(k; t) in Eq. (12.221) can be written in terms of this 
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particular solution as S(k;t) = S*(k;u(t)), with u(t) given by Eq. (I2.14p . This means 
that a sequence of static structure factors S*{k;u n ) generated by this simple exponential 
interpolating formula when the time t is given a sequence of values u n , say u n = nAu (with 
n — 0,1, 2, ...), will be identical to the sequence S(k; t n ) generated when the exact solution in 
Eq. (I2.22p is evaluated at a different sequence t n (n = 0, 1, 2, ...), provided that the times u n 
and t n are related by u n = J * n b(t')dt'. This observation greatly simplifies the mathematical 
analysis and the numerical method of solution of the full self-consistent theory under the 
particular conditions considered here. 

The solution / y(t) of Eq. (I2.2ip provides a dynamic order parameter in the sense that 
when it is infinite we can say that at that waiting time t the system remains ergodic, whereas 
if it is finite, we say that the system became dynamically arrested. A practical manner to 
use this criterion is to first construct a sequence of static structure factors S*{k;u n ) using 
Eq. (I2.23P for the uniform sequence u n = nAu (with n — 0,1, 2, ...). Each member of this 
sequence is then employed as the static input to solve self-consistently Eqs. (l2.16p -f l27L8|) . 
thus evaluating, using Eq. (I2.12p . a mobility sequence b{u n ). Since the sequence S*{k; u n ) is 
identical to the sequence S(k; t n ) provided that u n = J^ n b(t')dt', the mobility b{u n ) must be 
identical to b(t n ), and the corresponding time-sequence t n can be determined by means of 
the approximate recursive relation t n+ \ = t n + (Au)/b(t n ). If the dynamic arrest condition 
occurs along this process, i.e., if a value u^ a ' exists such that "f(u) (determined using S*(k; u) 
in Eq. (I2.2ip ) is infinite for u < u^ a ' and finite for u > u^ a \ then b(u) — > when u — > u^ a ' 
from below, and it is then not difficult to realize that the corresponding dynamic arrest 
time v- a > will diverge and = J °° b(t')dt' . The following numerical results illustrate the 
physical implications of this singular behavior. 

III. ILLUSTRATIVE APPLICATION. 

Let us now apply the theory just presented, to a concrete model system, namely, a 
dispersion of colloidal particles interacting through the hard-sphere plus attractive Yukawa 
pair potential expressed, in units of the thermal energy k B T = as 




(3.1) 
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FIG. 1: State space (4>,T*) of the hard-sphere plus attractive Yukawa model system (z = 20). The 
dotted line is the spinodal curve and the solid line is the dynamic arrest line calculated using Eq. 
(|2,2ip within the mean spherical approximation (MSA) for the equilibrium static structure factor 
S eq (k; (ft, T*). We consider an instantaneous quench process at t = from the ergodic initial state 
/ to the final state F near but slightly above the attractive glass "branch" of the dynamic arrest 
line. We also consider a second process, now to the point F' below this dynamic arrest line but 
still above the spinodal curve. 




The state space of this system is spanned by the volume fraction = ifna^j S /§ and the 
reduced temperature T* = K^ 1 , as illustrated in Fig. [TJ The equilibrium phase diagram of 
this system includes the gas and liquid disordered phases and crystalline solid phases. Here 
we will describe the equilibrium static structure factor S eq (k; <fr,T*) = \nS^ q (k; n, T*)]^ 1 of 
the disordered phases within the mean spherical approximation (MSA) 43[. Using this 
approximation and the compressibility equation 44 1 one can determine the spinodal curve 
of the gas-liquid transition by means of the condition 1/S eq (k = 0; <p,T*) = 0; the result is 
plotted in Fig. [T] for z = 20. 

Using the same MSA equilibrium static structure factor S eq (k; <fr,T*) in the equilibrium 
version of Eq. ( 12 . 2 1 [) we can scan the state space (<p,T*) to determine ^i eq at any point 
((f), T*) 30|. In this manner one locates the dynamic arrest transition line indicated by the 
solid curve of Fig. [TJ The region to the right and below this curve is thus predicted to 
correspond to dynamically arrested states. This figure focusses on the high-density glass- 
fluid-glass reentrance region that was first discovered using mode coupling theory j42|. We 
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now follow the approach introduced by Fofh et al. [14J in a simulation experiment on a very 
similar model system (a hard-sphere plus square- well fluid). Such experiment corresponds 
to suddenly quenching the system under isochoric conditions from an initial state (4>q,Tj) 
located in the fluid pocket of the reentrance (point / in Fig. [TJ, to a final state near the 
fluid- "attractive glass" transition line (either point F or point F' in Fig[T]). In the first case 
the end state (<po,Tp) lies slightly above the transition line, whereas in the second, the end 
state (<f>0jTp f ) lies in the region of arrested states. 

For this process we solve the general self-consistent system of equations in Eqs. f !2. 12[) - 
f !2.18p . The specific calculations are performed along the isochore 0o = 0.555 with initial 
temperature T/ = 0.159 and final temperature Tp = 0.0604. Fig. [2]illustrates the irreversible 
evolution of the static structure factor S(k; t w ) as a sequence of snapshots corresponding to 
five intermediate waiting times t w (from now on denoted by t w , rather than simply by t). We 
observe that the structure, initially described by S eq (k; 0o,T/), relaxes to the expected final 
value S eq (k;<f)Q,Tp), and that this process is faster at large wave-vectors, where it involves 
the appearance of stronger oscillations with k and a general shift of the maxima of S(k; t) 
to larger wave-vectors. To a large extent these features can be understood in terms even of 
the simple interpolating expression in Eq. f |2.23|) . 

The corresponding adjustment of the main peak of S(k; t w ) from its initial value 
S eq (k max ;(j)o,Tj) to its final value S eq (k max ; <j> , Tp) (< S eq (k max ; <f) ,Tj)) occurs, however, 
notoriously more slowly than at large wave-vectors and in an apparently non-monotonic 
manner, as illustrated in the inset of Fig. |2} which zooms on the evolution of the main peak. 
As observed there, as the system evolves, the maximum of S(k; t w ) moves to the right while 
decreasing in height to a value smaller than S eq (k max ; 4> ,T F ), bouncing back at later times 
to reach this final value. The origin of the predicted non-monotonic behavior can also be 
understood on the basis of the simple interpolation expression in Eq. (12.23p . which implies 
that S(k; t) will not change with waiting time for the wave-vectors k* at which the initial 
and the final static structure factors are already identical, S eq (k*;<po,Tj) = S eq (k*;<f)o,Tp). 
It is then not difficult to see that if the condition kinhc < k* < kmlx occurs, as it happens in 
our example, we shall observe this non-monotonic effect. 

A more interesting effect, which is perceptible in Fig. [2], but which is illustrated in 
more detail in Fig [3l is the evolution of S(k; t w ) at smaller wave-vectors. This refers to 
the emergence of a non-equilibrium low- A; peak that indicates the appearance of spatial 
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FIG. 2: Non-equilibrium evolution of the static structure factor S(k,t w ). The system, ini- 
tially equilibrated at (0o,2j) = (0.555,0.159), with S(k; t = 0) = S e i(k; <j> ,Tf) ((black) dot- 
ted curve), is instantaneously quenched at t w = to the final point (cf>Q,T F ) = (0.555,0.0604). 
The static structure factor then evolves continuously along a sequence of non-equilibrium 
values ((red) point-dashed lines) illustrated by the snapshots corresponding to t w /to = 
3.2, 60.87, 174.29, 945.39, 2023.54, and 4858.84, with t = [a 2 /D ]. Since the point F 
lies outside the dynamic arrest region, S(k,t w ) eventually attains its final equilibrium value 
S(k;t w = oo) = S eq (k; 4>o,T F ) ((blue) dashed curve). The main figure shows the resulting re- 
laxation process in a wide fc-range and the inset zooms on the evolution of the main peak of 



heterogeneities of average size \i(t w ) ~ 2n/ki(t w ), with k\ being the position of this emerging 
low- A; maximum. Fig [3ja provides a zoom on this effect in the case of the slightly deeper 
quench, now to the final state point F' in Fig. [TJ with temperature T F , = 0.0588 slightly 
below the dynamic arrest line. These heterogeneities may be associated with the appearance 
of voids whose average size and importance increase with waiting time, as suggested by the 
increasing height of the peak and by its shift to smaller wave- vectors observed as the system 
evolves. The emergence of this peak is associated with the vicinity of the gas-liquid spinodal 
region. In fact, it has the same origin as the low-/c peak that characterizes the process 



451 ]. even though in our case the final state (cj) ,T F ,) lies 



of early spinodal decomposition 
outside the spinodal region. 

As said above, this phenomenon is already observed in the shallower quench of Fig. |2j In 
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FIG. 3: Non-equilibrium evolution of the static structure factor S(k,t w ) for the deeper quench 
to the final state point F' . The system, initially equilibrated at (4>q,Tj) = (0.555,0.159), with 
S(k;t = 0) = S eg (k; 4>o,Tj) ((black) dotted curve), is instantaneously quenched at t w = to the 
final point(^, Tpi) = (0.555,0.0588) inside the dynamic arrest region. The static structure factor 
then evolves continuously along a sequence of non-equilibrium values ((red) point-dashed lines) 
illustrated by the snapshots corresponding to t w /to = 0.0, 1.23, 3.84, 140, 490, 975, 1265, and 
1522 tw /to)- Panel (a) focuses on the low-A; peak of S(k, t w ), and its inset shows the dependence 
of the position ki(t w ) of this low-fc peak on the waiting time t w (empty circles), with the solid line 
being the fit of the last few points with ki(t w ) rj (t w )~ a and a « ^. Panel (b) shows the behavior 
in a larger fc-regime, similar to Fig. [2j with its inset zooming on the evolution of the main peak. 

that case, however, although the system slows down considerably, the final structure of the 
irreversible evolution of S(k; t w ) is still the expected final equilibrium static structure factor 
S eq (k;(f) ,T F ), i.e., lim^^oo S(k; t w ) = S eq (k; (j) ,T F ) and the position k\(t w ) of this \ow-k 
peak decreases indefinitely. In contrast with that scenario, in the deeper quench illustrated 
in Fig[3j the final structure of the system is no longer S eq (k; </>o, Tpi); instead, the asymptotic 
\ong-t w limit of S(k; t w ) is given by S^ a '(k) = S*(k; u^ a '), where is the value of u at which 
the dynamic arrest condition is satisfied. This value is determined using the structure factor 
S*(k; u) of Eq. (I2.23j) as the structural input in Eq. (12. 21ft . as discussed at the end of Sec. 
UTI In Fig. |3ja we can compare the non-equilibrium arrested structure factor S^ a '(k) with 
the equilibrium structure S eq (k; <po,Tp) that would have been attained if no dynamic arrest 
condition had appeared along the equilibration process of S(k; t w ). 

In the same figure we also illustrate the evolution of S(k; t w ) towards its asymptotic limit 
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S( a \k) with a series of snapshots corresponding to a set of increasing waiting times. The 
most interesting feature revealed by these snapshots is the existence of an early evolution 
regime, in which S(k;t w ) evolves rather quickly towards the close neig hborhood of S^(k). 
As illustrated by these snapshots, this occurs within a finite waiting time d c) « 1500t . This 
early regime is followed by an asymptotic long-t w regime, in which the evolution of S(k;t w ) 
to actually reach the exact asymptotic value S^ a \k) becomes extremely slow and completely 
imperceptible in the scale of the figure. 

This is illustrated in the inset of Fig. |3^a), where we plot the evolution of the position 
ki(t w ) of the low-fc peak of S(k; t w ) for various waiting times between the last two snapshots 
of the main figure (i.e., 1265to < t w < 1522t ) • We notice that in this regime the last 
few data for ki{t w ) may be fitted approximately by a power law k\{t w ) ~ 10.22 x (t w )~ l l b . 

(c) 

In fact, the crossover waiting time t w can be estimated more accurately by the condition 
10.22 x (t^) -1 / 5 = k[ , with k^ = 2.34 being the asymptotic value of k\(t w ) corresponding 
to S {a) (k). This yields t$ « 1589t - The slow evolution regime t w > tw\ corresponding 
to asymptotically long times, cannot be observed, by definition, in the structural evolution 
illustrated in Fig. [3j It can, however, be observed in the evolution of the dynamic properties, 
as we discuss below. 

Let us emphasize the difference between the two quench processes just discussed (i.e., 
those involving the final state F or F'). For this, Fig. [3jb plots the evolution of S(k; t w ) for 
the latter process in the same manner as Fig. [2] does for the former. Let us point out that the 



quench simulated by Foffi et al. 



14j corresponds to the conditions illustrated in Fig. [3], i.e., 



to the process ending in the state F' just below the dynamic arrest line. We recall that in 
the process illustrated in Fig. [2] nothing prevents the evolution of S(k; t w ) from reaching the 
final structure factor S eq (k; <po, Tp), and this leads to the upturn of the peak illustrated in 
the inset of that figure. In contrast, as observed in the inset of Fig. [3jb, the main difference 
is that now the main peak of S(k; t w ) decreases but seems to stop evolving when t w reaches 
tw\ and this happens to occur before the upturn of the peak towards S eq (k; o ,Tp/) has a 
chance to develop. This is in agreement with what is observed in the simulated quench of 
Foffi et al., in which the main peak only decreases without exhibiting any upturn. On the 
other hand, our results in the inset of Fig. Ob also predict that the peak shifts slightly to 
the right, but in the simulation results such a shift is not appreciable. 

Similarly, in the report of the simulated quench of Foffi et al. lij no reference is made 
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FIG. 4: Theoretical predictions for the dependence of the intermediate scattering function 
F(k,r;t w ) on correlation time r for the quench to the final state F' corresponding to the waiting 
times t w /t = 0.0, 1.23, 3.84, 140, 490, 975, 1265, and 1522 (» t$/t ). In inset (a) and in the 
main figure the correlation time is scaled, respectively, by to an d by the relaxation time r a of the 
stretched exponential fit of the final relaxation of F(k,T;t w ). Inset (b) plots r a as a function of 
waiting time, with the line being the fit r a « i^ 7 (both times in units of to). 

to the low-fc peak predicted by our theory according to the illustrative results in Fig 0a. 
Thus, at this stage we cannot make a definitive statement on the level of a fine quantitative 
comparison between our theoretical predictions and the simulation results for the evolution 
of S(k;t w ), partially because of the differences in the model and in the conditions (volume 
fraction, for example) in which the quench was performed. While it is clearly desirable to 
carry out a systematic comparison on identical conditions, the agreement with important 
features observed in the simulation experiments is encouraging. 

Let us conclude this exercise by showing the irreversible evolution of the r-dependence 
of the intermediate scattering function F(k, r; t w ) for the quenching process I — > F' . This 
is presented in inset (a) of Fig. HI where the correlator f(k,r;t w ) = F(k,r;t w )/S(k;t w ) is 
plotted as a function of the correlation time r at representative waiting times corresponding 
to the snapshots of S(k; t w ) of Fig. [3J namely, t w /t = 0.0, 1.23, 3.84, 140, 490, 975, 1265, 
and 1522 (« tw /to). These results illustrate the fact that the decay of the temporal corre- 
lation of the fluctuations slows down notoriously as the system ages, developing a two-step 
relaxation: the initial /3-relaxation to an increasingly better denned plateau, followed by the 
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a-relaxation from this plateau. This is a typical behavior observed in the simulation and 
experimental studies of aging jl|, [§], liol-lol]. Another feature associated with aging is the 
superposition of the alpha relaxation at different waiting times on a single master curve, 
well-fitted by a stretched exponential function f(k,r;t w ) « A(k; t w ) exp[— (t/tq,)^]. Our 
theoretical results also exhibit this scaling property, as demonstrated in the main panel of 
Fig. HJ The exponent (3 is independent of t w (although it may depend on k). For the case 
illustrated in the figure we find (3 ~ 0.9. The a-relaxation time r a does depend on k and 
on t w , and the values of r a corresponding to each waiting time t w are plotted in the inset 
(b) of the same figure. At short times, these values are well fitted by a power law r a ~ tt 
characterized by the exponent z ~ 0.7. In the simulation experiment of Foffi et al. 14 1 
this scaling of the correlator is not fully apparent, although "in a crude tentative of data 
scaling", the authors report an exponent z ~ 0.38. At this point we should mention that, 
beyond detailed quantitative issues, the general predicted scenario illustrated in Fig. @] is 
completely similar to that reported in the simulated experiment of Foffi et al. 14J, which, 
in its turn, was found to be similar to that observed experimentally by Pham et al. js|. 

Regarding the low-fc peak predicted by our theory (see Fig[3]), let us notice that, although 
the final temperature of the quench is still above the spinodal temperature for this isochore, 
the asymptotic approach of S(k;t w ) to the non-equilibrium structure S^(k) is strongly 
suggestive of some form of arrested spinodal decomposition. In fact, preliminary calculations 
using our theory indicate that the scenario described in the main panel and the inset of 
our Fig. [3ja above, regarding this low-k peak and the phenomenon of arrested spinodal 
decomposition, is also predicted to occur at lower concentrations, in qualitative agreement 
with experimental observations (see Fig. 4.b and 4.c of Lu et al. pjj). Further comparisons 
and analysis lie, however, outside the scope of this illustrative presentation of the possible 
applications of the non-equilibrium SCGLE theory to the description of dynamic arrest 
phenomena, including aging, in instantaneously quenched uniform systems. 



IV. CONCLUDING REMARKS. 

In this manner, in section [TTT1 we have illustrated with a number of quantitative predictions 
for a specific model system (involving hard sphere plus short-ranged attractive interactions) 
the predictive nature of a generic theory of the non-equilibrium irreversible evolution of the 
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state of a homogeneous system subjected to a homogeneous and instantaneous quenching 
process. This theory is summarized by the self-consistent system of equations in Eqs. f !2. 1 2[) - 
(I2.18p . The time-evolving state of the system was described in terms of the static structure 
factor S(k; t w ) and of the r-dependence of the intermediate scattering function F(k, r; t w ) 
as a function of the waiting time t w after the quench. 

The specific process discussed corresponds to the sudden isochoric quench from an ini- 
tial fluid state (4>o,Tj) to a final state near the "attractive glass" transition. We observed 
that if the final state is also ergodic, the structure relaxes to its value equilibrium value 
S eq (k; 0o, Tp), whereas if the final state is in the dynamically arrested state, the struc- 
ture saturates asymptotically to a non-equilibrium value (k; (p ,T F ). In the latter case, 
S(k; t w ) develops a non-equilibrium low-fc peak that indicates the appearance of spatial het- 
erogeneities of average size \i(t w ) w 2ir/ki(t w ), with k\ being the position of this emerging 
low-A; maximum. The emergence of this peak is associated with the vicinity of the gas-liquid 
spinodal region. Regarding the evolution of the dynamics with aging time, the theory pre- 
dicts that the intermediate scattering function F(k, r; t w ) develops a two-step relaxations as 
the system ages. The theory also predicts the superposition of the alpha relaxation at differ- 
ent waiting times on a single master curve, well-fitted by a stretched exponential function, 
as observed in the simulation and experimental studies of aging. 

Let us stress that the theory proposed in section [Til however, is not limited to instanta- 
neous quench processes; in principle it is easily extendable to other quench "programs" by 
going one step back and use Eq. (12. lip instead of Eq. (12.131) . In this manner, a number 
of relevant questions could readily be addressed, such as the dependence of the aging of 
S(k; t w ) and F(k, r; t w ) on the quench protocol. Furthermore, in reality the theory of irre- 



versible relaxation in colloidal dispersions developed in Ref. [21], and summarized in Sec. 
[TTJ is not even limited to spatially homogeneous non-equilibrium states. The present work, 
however, was meant to provide the first exploratory application of this general theory in 
the simplest possible conditions. The specific results reported here suggest that this theory 
provides a qualitatively and quantitatively sound basis for the first-principles theoretical dis- 
cussion of the complex non-equilibrium phenomena associated with the aging of structural 
glass-forming colloidal systems. 
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